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INTRODUCTION 


In  the  US,  breast  cancer  is  the  second  m  ost  common  cancer  am  ongst  women  as  well  as 
the  second  leading  cause  of  death.  An  estim  ated  40,460  breast  cancer  deaths  occurred  in  2007.  1 
Mammographic  screening  can  help  detect  cancer  before  it  shows  sym  ptoms  and  dram  atically 
increase  survival  rates.  However,  the  f  alse  negative  rate  of  m  ammographic  procedures  is  in  10- 

30  %  rang  e,  thus  lim  iting  the  effectiveness  of  mammography  in  reducin  g  mortality."’  '  At  the 
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same  time,  the  false  positive  rate  is  estimated  to  be  11.5  %. 

In  mammography,  a  3D  volum  e  of  anatom  ical  structures  is  collapsed  into  a  2D  im  age 
plane.  The  resultant  im  age,  therefore,  is  a  consequence  of  projection  of  overlapping  anatomical 
structures  into  a  com  plex  2D  i  mage.  As  a  resu  It,  any  abnorm  ality  in  th  e  form  of  breast  cancer 
may  easily  get  hidden  behind  projections  of  a  norm  al  tissue  structures,  resulting  in  low 
sensitivity.  At  the  sam  e  time,  the  overlapping  st  ructures  may  also  mimic  the  pathology  that  the 
radiologists  are  looking  for,  leading  to  high  rate  of  false  positives.  An  imaging  techn  ique  which 
may  alleviate  the  lim  iting  factor  of  overlapping  anatomical  structures  and  at  the  same  time  take 
advantage  of  standard  m  ammography  im  aging  technique,  m  ay  ther  efore  prove  to  be  highly 
effective  in  breast  cancer  screening. 

The  Proposed  System:  Multi-plane  Correlation  Imaging  (MCI) 

MCI  is  a  realization  of  multi -projection  imaging  scheme  in  which  a  p  lurality  of  digital 
radiographic  images  of  the  sam  e  patient  is  acquire  d  within  a  short  inte  rval  of  time  from  slightly 
different  angles.  The  acquisition  protocol  is  similar  to  tomosynthesis  technique  with  a  difference 
that  th  e  im  ages  are  directly  analyzed  ins  tead  of  reconstructing  them  ,  thereby  avoiding 
reconstruction  artifacts.  These  images  ar  e  sim  ilar  to  projectio  n  im  ages  acqu  ired  in  s  tandard 
projection  technique,  except  that  each  of  the  a  ngular  projections  is  acquired  with  a  lower  dose 
level  such  that  the  total  patien  t  do  se  is  w  ithin  the  bounds  of  an  optim  ized  acquisition.  The 
acquired  images  are  then  processed  by  a  com  binatorial  technique  that  exploits  the  differences  in 
geometrical  perspec  tives  tha  t  dif  ferent  pro  jections  of  fer  to  e  stablish  sp  atial  cor  relation 
information  between  d  ifferent  angu  lar  p  rojections.  This  inf  ormation  is  used  to  id  entify  and 
positively  r  einforce  the  lesion  sig  nals  be  tween  dif  ferent  projec  tions,  thus  m  inimizing  the 
fundamental  lim  iting  factor  im  posed  by  anatom  ical  nois  e  on  the  d  etection  of  lesions.  A 
schematic  of  MCI  and  its  relationship  to  tomosynthesis  is  shown  in  Fig.l. 
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Multi-projection  Imaging  Scheme 

Multiple  Projection  Images 

Reconstruction  No  Reconstruction 


1  mage  slices 

Projection  images 

analyzed 

directly  analyzed 

Tomosynthesis 

/ 

MCI 

Fig.  1:  Schematic  of  Multi-plane  Correlation  Imaging  (MCI)  illustrating  its  lineage. 


Clinical  Significance  of  MCI 

MCI  m  eets  the  unm  et  need  in  diagnostic  radiology  for  a  technique  that  enables 
radiologists  to  faithfull  y  differentiate  cancer  from  nor  mal  tissue  and  sh  ow  ca  ncers  hidden  by 
superimposition  of  anatom  ical  structure.  W  hen  combined  with  a  rapid  acquisition  digital  x-ray 
detector,  MCI  m  ay  be  im  plemented  on  m  odified  radiographic  im  aging  system  s  in  a  cost- 
effective  way,  and  incorporated  to  yield  hi  gh  patient  throughput.  P  ractically,  MCI  can  be 
implemented  in  different  form  s  including  scrolli  ng  the  im  ages  m  anually  or  in  cine  m  ode, 
stereoscopic  display  of  projections  images,  or  computer-aided  analysis  of  the  images.  Overall,  if 
perfected  and  im  plemented  using  optim  ized  da  ta  acquisition  schem  e,  MCI  could  provide  a 
simple  and  accurate  can  cer  diagnosis  over  projec  tion  imaging  techniques  at  lower  patien  t  dose 
levels. 

Compared  to  tom  osynthesis,  MCI  technique  could  provide  a  sim  pier,  fast  er  and  1  ower 
cost  alternative  with  comparable  or  better  performance  for  the  detection  of  cancer,  potentially  at 
lower  dose  levels.  Furtherm  ore,  the  absence  of  the  need  to  reconstruct  im  ages,  and  thus  the 
associated  artifacts  inherent  in  tomosynthesis  imaging,  provides  the  MCI  technique  an  immediate 
advantage  in  terms  of  the  diagnostic  image  quality.  Moreover,  as  compared  to  tomosynthesis  in 
which  as  many  as  50-80  slices  may  need  to  be  reviewed  depending  on  the  size  of  the  patient,  the 
MCI  technique  m  ay  entail  inves  tigating  a  significantly  fewer  num  her  of  im  ages,  resulting  in 
higher  throughput  and  potentially  improved  confidence  or  even  accuracy  of  a  radiologist’s 
decision. 
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We  also  envision  the  use  of  MCI  as  an  adj  unct  to  tomosynthesis  imaging,  with  no  added 
patient  dose.  The  CAD  technique  developed  for  MCI  m  ay  be  applied  to  the  projection  im  ages 
acquired  in  the  process  of  tomosynthesis  im  aging,  providing  the  radiologist  with  a  second 
diagnostic  opinion.  Thus,  as  with  tomosynthesis,  besides  mammography  and  chest  imaging,  MCI 
could  be  used  in  a  wide  ranging  applicati  ons  including,  but  not  lim  ited  to,  brachytherapy, 9  hand 
arthritis,10  dental,11  and  radiotherapy  applications.12 

Optimization  of  MCI 

In  MCI  several  possible  com  binations  of  its  data  acquisition  param  eters  such  as  the 
acquisition  dose  lev  el,  the  num  her  of  angular  acqui  sitions,  and  the  total  a  ngular  sp  an  of  those 
acquisitions  m  ay  be  used.  Not  a  11  of  thes  e  com  binations,  however,  are  optim  al  -  c  ertain 
combination  of  these  parameters  may  quickly  distort  pathological  indicators,  making  it  important 
to  derive  an  optimum  acquisition  scheme  that  would  make  MCI  maximally  effective. 

While  some  initial  studies  on  optimization  of  tomosynthesi  s  have  been  conducted  they 
cannot  be  directly  applied  to  MCI.  This  is  because  most  of  the  current  optimization  protocols  are 
not  based  on  m  aximizing  the  quality  of  output  images  -  which  is  an  im  portant  f  actor  in 
characterizing  the  clinical  utility  of  a  system  .  The  optimization,  therefore,  should  focus  on  the 
dependence  of  i  maging  system  perfonn  ance  on  system  configuration  and  derive  a  specific 
configuration  that  delivers  best  system  performance.  An  optimized  data  acquisition  scheme  will 
yield  superior  diagnostic  accuracy  of  MCI,  particularly  importa  nt  for  the  early  detection  of 
cancer.  Furthermore,  optimizing  acquisition  parameters  will  allow  to  fully  realize  the  potential  of 
MCI  in  delivering  superior  diagnostic  perform  ance  at  a  dose  level  lower  than  that  delivered  in 
standard  procedures.  This  would  be  in  line  with  the  “as  low  as  reasonably  achievable”  (ALARA) 
dose  principle  governing  cancer  screening. 


Proposed  Work 

In  this  work,  we  inves  tigated  the  f  easibility  of  a  new  optim  ized  image  acquisition  and 
processing  approach,  n  amely.  Multi-projection  Correlation  Imaging  (MCI).  The  hypothesis  of 
this  study  was  that  the  detection  of  cancer  can  be  im  proved  by  m  inimizing  the  undesirable 
influence  of  overlying  anatomical  noise  by  ham  essing  correlation  infonnation  between  multiple 
projections  of  the  sa  me  patient  in  an  optim  ized  fashion.  Grounded  on  this  hypothesis,  this 
research  aim  ed  to  (a)  theoretica  lly  dem  onstrate  the  clin  ical  advan  tage  in  using  MCI  over 
standard  projection  technique  as  a  new  dia  gnostic  tool  for  im  proved  cancer  detection,  (b) 
develop  and  validate  a  framework  to  optimize  the  image  acquisition  process  in  MCI,  (c)  use  this 
framework  to  establish  an  optim  ization  rule  for  MCI  and  its  o  ther  multi -p  rojection  im  aging 
counterpart,  namely  tomosynthesis,  and  finally  (d)  compare  clinical  pe  rformance  of  MCI  with 
tomosynthesis  to  examine  its  clinical  utility. 
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BODY 


Specific  Aim  1:  Determine  the  set  of  acquisition  parameters  for  an  MCI  study.  (Months  1-9) 

Progress:  Work  for  this  task  h  as  been  accomplished  and  was  detailed  in  the  last  an  nual  report. 
The  goal  of  this  task  was  to  acquire  images  that  approximate  mammographic  backgrounds  which 
could  then  be  used  for  analysis  in  developing  m  ethodologies  for  the  other  specific  aim  s  of  this 
study. 

Conclusions  of  Task  1 


Sub-optimized  im  plementation  of  MCI  can  p  otentially  c  ompromise  it  s  ma  ximum 
achievable  diagnostic  perform  ance.  In  this  ta  sk,  we  developed  an  algorithm  ic  observer-based 
framework  to  assess  the  im  pact  of  various  acquisition  parameters  of  MCI  performance,  namely, 
dose,  the  number  of  projections,  and  the  angular  span.  The  study  dem  onstrated  the  interplay  of 
anatomical  and  quantum  noise  in  the  overall  perform  ance.  The  following  conclusions  were 
drawn  from  this  phase  of  the  study: 

1)  Increasing  the  number  of  projections  while  keeping  the  overall  dose  and  angular 
span  constant  decreased  the  perfonnance  of  MCI. 

2)  Increasing  the  angular  span  and  acquisi  tion  dose  level  improved  the  m  aximum 
obtainable  AUC. 

3)  The  number  of  projections  required  to  maximize  perfonnance  was  found  to  be 
linearly  related  to  the  angular  span.  This  num  her  was  found  to  be  independent  of  the  acquisition 
dose  level.  The  best  clinical  perfonn  ance  wa  s  obtained  when  the  angu  lar  sep  aration  between 
each  projection  was  approximately  2.75°. 

4)  Finally,  the  results  revealed  th  at  the  peak  performance  for  MCI  at  the  clinically 
relevant  dose  levels  of  one-  a  nd  two-view  m  ammography  was  achie  ved  at  15  -  17  projections 
spanning  an  angular  arc  of  ~45°,  the  widest  angle  tested  in  this  study. 

The  outcome  of  this  task  was  published  in  two  ar  tides  the  journal  of  Medical  Physics  in  2007 
and  2008: 

(a)  Chawla  A. ,  Samei  E.,  Saunders  R.,  Abbey  C.,  Del  ong  D.,  Effect  of  dose  reduction  on  the 
detection  of  mammographic  lesions:  A  mathematical  observer  model  analysis,  Medical  Physics 
34:3385-3398,2007. 

(b)  Cha  wla  A. ,  Sa  mei  E.,  Saunders  R.,  Lo  J.,  Baker  J  .,  A  m  athematical  m  odel  platfonn  for 
optimizing  a  multi-projection  breast  imaging  system,  Medical  Physics  35:  1337-1345,  2008. 
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Specific  Aim  2:  Extend  single-view  CAD  processi  ng  m  ethods  used  for  conventional 
mammography  for  MCI  implementation  using  a  multi-plane  correlation  rule.  (Months  10-21) 

Progress:  Work  for  this  task  has  been  com  pleted  as  well  and  was  detailed  in  the  last  annual 
report.  A  CADe  processor  developed  earlier  for  standard  projection  technique  was  extended  to 
take  advantage  of  the  MCI  configuration.14 

Conclusions  of  Task  2 


A  new  CADe  processor  was  developed  for  multi-projection  Correlation  Imaging  (Cl)  that 
takes  advantage  of  the  g  eometrical  correlation  information  accrued  from  the  ava  ilable  multiple 
projections  to  im  prove  specificity  of  the  Cl  system.  The  perfor  mance  of  CADe  proc  essor  was 
found  to  be  robust  as  it  is  successfully  able  to  locate  the  suspected  lesions. 

The  work  resulting  from  this  task  was  subm  itted  to  the  journ  al  of  Academic  Radiology  and 
has  now  been  accepted  for  publication.  The  s  ubmitted  m  anuscript  h  as  been  attached  in  the 
appendix. 

-  Cha  wla  A.,  Samei  E.,  Lo  J.Y,  Baker  J.,  Toward  Optim  ized  Acquisition  Sche  me  f  or  Multi¬ 
projection  Correlation  Imaging  of  Breast  Cancer,  In  print,  Academic  Radiology,  2008. 


Specific  Aim  3:  Fine-tune  the  optimized  formalism  of  step  1  for  CAD  algorithm  and  evaluate 
its  performance. 


Progress:  Work  f  or  this  task  ha  s  been  completed.  The  optim  ization  results  o  btained  with 
observer  model  processor  in  step  1  were  com  pared  and  sub  stantiated  with  tho  se  obtained  with 
CADe  processor. 

Task  3.1:  Apply  the  CAD  algorithm  on  multi-projection  images  of  breast  specimens  to 

test  the  optimality  of  parameters  established  in  task  1.  (Months  22-25) 

To  meet  the  goals  of  this  aim,  the  acquisition  parameters  were  systematically  changed  and  the 
CADe-based  performance  measured  for  different  settings  of  those  parameters.  These  results  were 
compared  to  the  observer  model-based  performance  to  confirm  the  optimality  of  MCI  scheme. 

The  specific  methods  developed  to  accomplish  this  task  will  be  detailed  here. 
a.  Image  Database 

The  study  employed  a  database  of  image  sets  from  106  subjects  recruited  for  our  ongoing 
tomosynthesis  clinical  trial.15  Each  image  set  consisted  of  25  i  mages  of  a  single  breast  acquired 
about  the  CC  or  MLO  orientation  from  25  different  but  fixed  angular  positions  unifonnly  spaced 
in  steps  of  1.8  °  over  a  -45  °  arc.  A  prototype  clin  ical  m  ulti-projection  sys  tern,  Siem  ens’ 
Mammomat  Novation  TOMO  (Fig.  2)  was  used.  16  The  im  ages  were  acquired  at  kVps  ranging 
between  28  and  30,  with  a  total  glandular  dose  less  than  that  delivered  in  a  standard  two-view 
screening  p  rocedure.  A  11  cases  were  interpreted  by  one  of  fi  ve  dedicated  breast  im  aging 
radiologists  to  be  normal  (without  any  lesions). 
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Translating  x- 
ray  source 
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Fig.  2:  Schematic  of  acquisition  for  multi-projection  breast  Correlation  Imaging  (Cl). 

Front  view  (left);  side  view  (right) 

Next,  a  database  of  lesion  present  im  ages  was  generated.  Toward  that  end,  53  out  of  the 

available  106  cases  in  the  database  of  normal  clinical  images  were  supplemented  with  projections 

from  a  simulated  3D  3  mm  lesion.17  The  lesion  was  simulated  to  be  located  at  the  center  of  the  breast 

a  distance  of  3  c  m  above  the  detector  plane.  The  projections  of  this  lesi  on  on  the  detector  were 

simulated  for  all  the  25  different  tube  positi  ons  and  embedded  into  corresponding  angular 

projections  of  the  subject  images.  Thus  there  were  two  datasets  each  of  53  s  ubjects,  one  with  lesion 

absent  and  the  other  with  lesion  present.  The  contrast  of  the  lesion  was  set  assuming  a  heterogeneous 

breast  (50%  glandular/50%  adipos  e  tissue,  representin  g  an  av  erage  breast  com  position)  and 

accounting  for  the  acquisition  kVp,  target/filter  combination,  breast  thickness,  anode  type,  and 
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appropriate  scatter  fractions. 

b.  Optimization  of  Data  Acquisition  and  Comparison  with  Observer-model 
Results 

To  optimize  the  acquisition  scheme,  the  components  of  acquisition,  namely,  the  number  of 
projections  and  their  angular  span  were  sy  stematically  changed  within  2-25  and  3.6  0  -  44.8  0 

range,  respectively,  to  investigate  which  one  of  the  m  any  possible  com  binations  yielded  the 
highest  diagnostic  performance. 

Area  under  the  ROC  curve  (AUC)  was  derived  fir  om  the  datasets  w  ith  and  without  the 
embedded  lesion.  Each  case  in  the  tw  o  datasets  was  processed  with  the  C  ADe  processor  described 
earlier  to  yield  a  corresponding  2D  contour  m  ap.  The  likelihood  of  the  presen  ce  of  the  em  bedded 
lesion  was  examined  via  a  correlation  m  atching  of  the  expected  signal  with  the  signal-present  and 
with  the  signal-absent  2D  contour  m  ap.  The  value  obtained  by  this  signal-matching  step  served  as 
decision  variables  based  upon  whic  h  the  probability  distribution  functions  (  pdf)  of  the  signal- 
absent  and  signal-present  decisi  on  variables  were  computed.  Finally,  non-param etric  ROC  curves 
were  derived  by  thresholding  the  pdfs,  and  AUC  subsequently  co  mputed  by  the  trapezo  idal 
method. 
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AUC  was  com  puted  using  location  and  signal  known  exactly  (SKE)  paradigm  similar  to 
the  procedure  used  by  m  athematical  observe  r  m  odel  processors.  Th  erefore,  although  not 
equivalent,  the  AUC  values  obtained  by  the  tw  o  processors  are  com  parable.  Most  importantly, 
they  provide  a  single  platform  on  which  to  comp  are  the  p  erformance  results  of  o  bserver-model 
may  be  compared  with  those  from  CAD. 


Results  of  Task  3.1 


Fig.  3  show  s  the  variation  of  AUC  with  the  number  of  projections  spanning  different 
angular  arcs.  At  each  angular  range,  the  AUC  values  increase  with  the  increase  in  the  num  her  of 
angular  projections  and  then  appear  to  approach  an  asym  ptote.  The  num  her  of  projections  at 
which  the  A  UC  values  peak  depends  on  the  angul  ar  span.  The  highest  AUC  is  obtained  at  the 
maximum  angular  span  of  44.8°  with  7  projections. 

The  trend  in  the  variation  of  AUC  values  de  lineate  the  role  of  different  com  ponents  of 
acquisition  schem  e  in  the  final  diagnostic  perf  ormance  of  a  m  ulti -projection  im  aging  sys  tern. 
These  trends  indicate  that  the  optimum  number  of  projections  for  a  multi-p  rojection  imaging 
system  may  be  in  the  7-10  range  for  an  angular  span  of  44.8  °.  Most  noteworthy,  the  observer 
model  results  (reported  in  the  last  annual  report  and  reproduced  here  in  Fig.  4  for  convenience) 
show  a  sim  ilar  trend  in  perfor  mance  where  the  m  aximum  detec  tability  of  an  em  bedded  lesion 
was  found  to  maximize  with  between  10-17  angular  projections  for  an  angular  span  of  44.8°. 
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AUC  val  ues  AUC  val  ues  AUC  val  ues 


Number  of  Projections  Number  of  Projections 


Fig.  3:  Area  under  ROC  curves  as  a  function  of  the  number  of  projections  spanning 
different  angular  ranges  (specific  values  shown  in  the  legends)  in  a  Multi-plane  Correlation 
Imaging  setup.  AUCs  indicate  the  detectability  of  a  simulated  mass  embedded  into  each 

projection. 
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Fig.  4:  Variation  of  AUC  for  different  number  of  angular  projections  spanning  for  representative 
angular  spans  in  the  7.5°-44.8°  range  using  a  mathematical  observer  model.  These  results 
confirm  the  optimization  results  obtained  from  the  CADe  processor  (shown  in  Fig.  3). 


Task  3.2:  Compare  the  detection  performance  of  extended  CAD  algorithm  for  MCI 
with  that  of  an  available  single-view  CAD  algorithm  used  in  conventional 
mammography.  (Months  26-28) 

This  task  w  as  accomplished  as  a  part  of  task  3  .1  in  which  the  performance  of  CADe  for 
MCI  was  determ  ined  at  different  acquisition  settings,  in  eluding  that  for  singl  e-view  conventional 
mammography.  Specifically,  the  di  agnostic  perform  ance  was  m  easured  in  tenn  s  of  two 
performance  indices,  first  as  the  ratio:  True  Positives/(True  Positives  +  False  Positives).  This  ratio, 
tenned  Positive  Predictive  Index  (P  PI),  is  a  m  easure  of  the  true  positive  lo  cations  as  a  fraction  of 
the  to  tal  nu  mber  of  ide  ntified  loca  tions  pe  r  im  age  s  et  a  nd  i  s  e  asily  de  rivable  fir  om  t  he  FRO  C 
curves.  These  values  were  then  averaged  across  ah  the  cases  for  each  possible  com  bination  of  the 
number  of  projections  and  angular  range. 

Fig.  5  shows  the  variation  in  the  average  posi  five  predictive  index  (PPI)  with  the  number 
of  projections  within  6  angular  spans  in  the  7.5  °-44.8°  range.  Most  noteworthy,  at  each  angular 
range,  the  P  PI  values  first  increas  es  with  increase  in  the  number  of  projections.  This  suggests 
that  CADe  for  MCI  outperform  ed  CADe  for  single-view  conventional  mammography.  The  PPI 
values  furth  er  increas  e  with  increas  e  in  the  n  umber  of  projections  and  then  decrease  with  an 
increase  in  the  num  ber  of  projections,  m  aximizing  at  a  value  that  is  dependent  on  the  angular 
span.  The  maximum  PPI  is  obtained  for  10  projections  spanning  an  angular  arc  of  44.8°. 
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Fig.  5:  Average  Positive  Predictive  index  [TP/(TP  +  FP)]  as  a  function  of  the  number  of 
projections  spanning  different  angular  ranges  (specific  values  shown  in  the  legends)  in  a  multi¬ 
projection  Correlation  Imaging  setup.  TP~True  Positive  findings;  FP~  False  Positive  findings  per 

patient  case. 


Conclusions  of  Task  3 


In  conclusion,  the  perfor  mance  of  the  CADe  system  was  com  puted  at  different  data 
acquisition  settings  towards  optimizing  the  geometry  of  image  acquisition  of  MCI.  Both  the  CADe 
and  observer  m  odel  results  (reported  earlier)  show  a  mutually  re  inforcing  tr  ends  in  th  e 
performance  of  Cl  as  a  function  of  the  differe  nt  acquisition  com  ponents,  and  confirm  that  the 
maximum  performance  may  be  obtained  with  7-17  projections  for  an  angular  span  of  ~45  °.  Most 
noteworthy,  the  optimization  framework  reported  here  is  based  on  two  se  parate  clinically-relevant 
processors  and  thus  could  be  used  to  potentially  improve  the  clinical  efficiency  of  MCI. 
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Specific  Aim  4:  Evaluate  the  clinical  performance  of  MCI  and  its  improvement  over  the 
conventional  single-view  CAD  system  in  the  detection  of  breast  cancer  using  a  limited  number 
of  human  subject  cases. 

Work  for  this  task  was  substituted  with  a  m  ore  clinically  relevant  task  of  comparing  MCI 
with  breast  tom  osynthesis.  Th  is  is  becaus  e  breast  tom  osynthesis  is  a  new  b  reast  im  aging 
technique  that  has  show  n  promising  results  in  imp  roving  accuracy  of  breas  t  cancer  detection  in 
recent  clinical  trials  and  may  possibly  receive  PMA  FDA  approval  for  clinical  use  in  2009.  Thus, 
to  evaluate  clinical  utility  of  MCI  for  eventual  clinical  adoption,  it  was  paramount  to  compare  its 
performance  to  breast  tomosynthesis. 

Towards  meeting  the  goal  of  this  task,  the  framework  based  on  projection  im  ages  reported 
earlier,  was  extended  to  include  to  mosynthesis  reconstruction.  The  objectives  of  this  task  were 
threefold:  (a)  to  investigate  th  e  role  of  acquistion  param  eters  on  the  diagnostic  quality  of  breast 
tomosynthesis  at  clinically  relevant  dose  levels,  (b)  to  derive  a  comprehensive  optimization  rule 
for  tomosynthesis  in  terms  of  a  specific  com  bination  of  the  acquis  ition  components  that  renders 
the  best  available  diagnostic  information,  and  finally,  (c)  to  quantitatively  compare  the  optimized 
performances  of  MCI  and  tomosynthesis. 

This  task  was  the  major  f  ocus  of  Pi’s  ef  forts  in  the  last  year,  and  has  now  been 
accomplished.  The  specific  methods  and  framework  developed  for  this  phase  to  meet  the  goals  of 
this  task  will  be  detailed  here. 

I.  Materials  and  Methods 
A.  Image  database 

Images  from  nice  m  astectomy  specim  ens  th  at  were  used  in  Task  1  to  analy  ze  th  e 
performance  of  MCI  were  em  ployed  to  investig  ate  the  perform  ance  of  tom  osynthesis.  To 
incorporate  the  specific  configuration  of  to  mosynthesis,  reconstructed  slice  volum  e  was 
generated  for  each  m  astectomy  specimen.  To  that  end,  the  projection  im  ages  with  and  without 
supplemented  lesions  w  ere  reconstructed  using  Siemens’  proprietary  filtered  back  projection 
(FBP)-based  reconstruction  algorithm,  Tom  oEngine.19  Variab  le  im  age  slice  vo  lumes  were 
generated  depending  on  the  compressed  size  of  the  specimens  to  generate  fixed  slice  thickness  of 
1  mm. 


B.  Evaluation  of  Tomosynthesis  Detection  Performance 

For  evaluating  tom  osynthesis,  the  in-focus  slice  corresponding  to  the  lesion  plane  3cm 
above  the  detector  surface  was  extracted  from  the  reconstructed  image  slices  for  each  specim  en. 
The  specific  param  eters  for  reconstruction  we  re  dictated  by  the  acquisition  configuration 
considered.  The  84  lesion  supplem  ented  ROIs/sp  ecimen  were  then  extracted  from  the  slice 
resulting  in  756  ROIs  for  all  the  specimens.  These  ROIs  were  analyzed  for  the  presence  of  signal 
using  the  mathematical  model-based  methodology  similar  to  that  used  for  Cl  .  At  the  end  of  this 
analysis,  a  single  ROC  curve  was  derived  corresponding  to  the  overall  detectability  of  lesions  in 
a  reconstructed  slice. 
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C.  Optimization  Framework 


The  optim  ization  fram  ework  was  sim  ilar  to  that  used  for  MCI,  except  that  it  was 
extended  to  include  reconstruc  tion  per  tom  osynthesis  conligur  ation.  F  ig.  6  provides  a  visual 
illustration  of  this  multi-factorial  optimization  scheme.  The  specific  combination  of  acquisition 
components  that  yielded  the  m  aximum  detect  ion  performance  was  deem  ed  the  optim  ized 
acquisition  parameters  set  for  tomosynthesis. 


Fig.  6:  Schematic  of  the  optimization  space  used  in  this  study  to  analyze  MCI. 


II.  Results 

Fig.  7  shows  variation  in  AUC  with  the  num  her  of  angular  projections  at  different  angular 
ranges  for  MCI  and  reconstructed  slices  (tom  osynthesis),  respectively,  under  iso-im  age  dose 
condition  with  dose  lev  el  of  each  projection  was  fixed  at  1/25  th  fraction  of  typical  single-view 
mammography.  (It  m  ay  be  noted  that  the  results  fo  r  MCI  have  been  report  ed  previously  in  the 
last  annual  report  but  have  been  reproduced  here  for  convenience  such  that  a  direct  com  parison 
of  MCI  and  tomosynthesis  may  be  drawn.) 

For  both  im  aging  modes,  regardless  of  the  angu  lar  span,  the  AUCs  f  irst  increased  with  an 
increase  in  the  number  of  projections  but  leveled  off  beyond  a  certain  number  of  projections.  The 
maximum  value  of  AUC,  however,  increased  with  an  increase  in  the  a  ngular  span.  Projection 
images  yielded  a  slightly  higher  m  aximum  AUC  than  reconstruc  ted  slic  es  indicating  slight 
benefit  in  the  detection  perform  ance  in  using  just  the  projection  im  ages  over  tom  osynthesis. 
Most  im  portantly,  the  peak  perform  ance  fo  r  both  im  aging  m  odes  was  between  15  and  20 
projections  for  an  angular  span  of  about  45°. 

Figs.  8  show  reconstructed  in- focus  slice  of  a  m  astectomy  specimen  under  iso-im  age  dose 
condition  of  single-  view  m  ammography  dose  level.  The  slice  was  generated  using  5  (a),  13  (b), 
and  25  (c)  projections  spanning  a  fixed  angular  range  of  44.8°.  There  is  a  notable  improvement  in 
the  detectability  of  les  ions  as  th  e  number  of  pr  ojections  increases  from  5  to  13;  however,  the 
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detectability  does  not  appear  to  improve  by  incr  easing  the  num  ber  of  projections  to  25.  This 
provides  a  visual  evidence  of  the  saturation  in  AUC  values  shown  in  Fig.  7b. 
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Fig.  7:  Variation  of  AUC  with  number  of  projections  for  Cl  (a)  and  tomosynthesis  (b),  at  different 
acquisition  dose  levels  under  iso-image  dose  condition  (dose  of  each  projection  was  fixed  at  l/25th  fraction 
of  typical  single- view  mammography).  The  angular  spans  of  these  projections  were  in  the  7.5°  -  45°  range. 
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Fig.  8:  In-focus  tomosynthesis  slice  at  the  central  plane  of  lesions  embedded  in  a  mastectomy  specimen.  The  slice  was 
reconstructed  using  5  (a),  13  (b),  and  25  (c)  projections  under  iso-image  dose  condition  (dose  of  each  projection  used 
for  reconstruction  was  fixed  at  1/25*  fraction  of  typical  single-view  mammography).  The  detectability  of  the  lesion 
evidently  improves  from  5  to  13  projections,  and  then  appears  to  remain  steady  beyond  13,  thus  visually  confirming 

saturation  in  AUC  values  shown  in  Fig.  7b. 
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Figs.  9  and  10  show  the  variation  in  AUC  with  the  angular  span  and  number  of  projection 
for  projection  images  and  reconstructed  slices  under  fixed  dose  le  vels  (iso-study  dose  condition) 
at  total  dose  levels  equivalent  to  that  of  half  and  full  single -view  mammography  dose  level, 

respectively.  Both  projection  images  and  tomosynthesis  showed  improvement  in  diagnosis  when 
information  from  multiple  images  was  combined,  confirming  the  benefit  of  both  methods  over 
standard  m  ammography.  However,  for  all  angula  r  spans,  the  AUC  fir  st  increased  and  then 
decreased  as  the  num  her  of  projections  was  incr  eased.  The  num  her  of  pr  ojections  at  which  the 
AUC  values  peak  was  dependent  on  the  angul  ar  span.  Most  noteworthy,  the  m  aximum  AUC 
value  for  both  projection  im  ages  a  nd  reconstructe  d  slices  was  obtained  at  an  angular  span  of 
44.8°  with  15  -  17  projections.  This  su  ggests  that  current  implementations  of  breast 

tomosynthesis  with  25  projections  may  be  sub-optim  al.  Reconstructed  slices  were  also  noted  to 
have  slightly  lower  perfor  mance  than  projectio  n  im  ages,  especia  lly  at  the  ha  If  single-v  iew 
mammography  dose  level  (shown  in  Fig.  9b),  possibly  due  to  reconstruction  artifacts  that  lim  it 
the  efficiency  of  tomosynthesis.  The  higher  performance  of  projection  images  is  thus  potentially 
indicative  of  the  maximum  achievable  perfonnance  via  tomosynthesis. 

Fig.  1 1  provides  a  visual  evidence  of  change  in  detectability  of  the  em  bedded  lesions  in 
the  central  slice  under  iso-study  dose  condition  at  dos  e  levels  equivalent  to  that  of  single-  and 
two-view  mammography,  respectivel  y.  The  number  of  angular  proj  ections  used  to  reconstruct 
the  slice  was  varied  from  5  (a)  to  13  (b)  to  25  (c)  spanning  a  fixed  angular  arc  of  44.8  °.  At  both 
the  dose  levels,  the  detectability  of  the  lesions  evidently  improved  from  5  to  13  projections  but 
decreased  when  the  number  of  projections  was  increased  to  25. 

Figs.  12  and  13  depict  that  finding  for  pr  ojection  im  ages  and  tom  osynthesis.  Fig.  12 
shows  the  num  her  of  pr  ojections  that  yield  m  aximum  AUC  at  differe  nt  angular  ranges,  while 
Fig.  13  shows  the  corresponding  AUCs  at  each  of  those  angular  ranges.  The  maximum  AUC  is 
obtained  using  a  44.8°  angular  span  and  15-17  projections.  As  44.8°  was  the  maximum  angular 
span  tes  ted  in  this  stud  y,  it  is  exp  ected  that  a  wider  ang  ular  sp  an  might  yield  even  higher 
performance.  The  slope  of  the  linear  fit  in  Fig.  12  reveals  that  for  projection  im  ages  and 
reconstructed  slices,  the  optim  urn  angular  separation  that  realizes  maxim  um  perfor  mance  is 
approximately  2.75°. 
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Fig.  9:  Variation  of  AUC  with  the  number  of  projections  for  projections  images  (a)  and  reconstructed  slices  (b)  under 
iso-study  dose  conditions  at  different  angular  ranges  in  7.5°  -  45°  range.  The  total  dose  level  was  fixed  to  half  that  of 

single- view  mammography. 
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Fig.  10:  Variation  of  AUC  with  the  number  of  projections  for  projection  images  (a)  and  reconstructed  slices  (b)  under 
iso-study  dose  conditions  at  different  angular  ranges  in  7.5°  -  45°  range.  The  total  dose  level  was  fixed  to  that  of 

single- view  mammography. 
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Fig.  11:  In-focus  tomosynthesis  slice  at  the  central  plane  of  lesions  embedded  in  a  mastectomy  specimen.  The  slice 
was  reconstructed  using  5  (a),  13  (b),  and  25  (c)  projections.  The  total  dose  was  fixed  at  that  of  single-view 
mammography  regardless  of  the  number  of  projections.  The  detectability  of  the  lesion  evidently  improves  from  5  to 
13  projections,  and  then  appears  to  deteriorate  at  25  projections  due  to  increased  noise  per  projection,  thus  visually 

confirming  the  roll  off  in  the  AUC  values  shown  in  Fig.  10b. 
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Projs  =  0. 38*Ang  Range 


Fig.  12:  The  number  of  projections  per  angular  range  that  yield  maximum  AUCs  for  projection  images  (a) 
and  reconstructed  slices  (b).  These  are  plotted  for  different  dose  levels  (denoted  n  the  legends  as  the 

multiples  of  that  of  single-view  mammography). 
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KEY  RESEARCH  ACCOMPLISHMENTS  IN  YEAR  2006-09 


In  th  is  research  work,  we  assessed  the  f  easibility  of  MCI  techn  ique,  investigated  its 
optimization  for  m  aximum  diagnostic  output,  and  fi  nally  demonstrated  its  utility  for  im  proved 
breast  cancer  detection. 

Three  of  th  e  four  specific  aim  s  outlined  in  this  study  were  successfully  m  et.  The  final 
goal  of  limited  MCI  clinical  trial  w  as  substituted  for  a  m  ore  important  task  of  comparing  th  e 
clinical  perf  onnance  of  MCI  with  that  of  tomo  synthesis.  Such  a  comparison  was  considered 
critical  since  it  would  have  a  bearing  on  establishing  clinical  utility  and  thus,  ultimately,  potential 
clinical  adoption  of  MCI. 

More  specifically,  in  this  st  udy,  a  theoretical  m  odel  was  first  developed  to  determ  ine  the 
diagnostic  in  for  mation  content  of  projection  im  ages  in  MCI  us  ing  a  m  athematical  observer. 
Using  this  model,  a  multi-f  actorial  task-based  framework  was  developed  to  optim  ize  the  image 
acquisition  of  MCI  using  existing  low-dose  clini  cal  data.  The  fra  mework  was  further  validated 
using  a  CADe  processor.  Performance  of  MC  I  was  evaluated  on  m  astectomy  specim  ens  at 
clinically  relevant  doses  and  further  compared  to  tomosynthesis. 

Application  of  the  theo  retical  m  odel  on  low-d  ose  clinical  data  s  howed  that  peak  MCI 
performance  m  ay  be  obtained  with  15-17  projec  tions.  CAD  results  confirm  ed  si  milar  trends. 
Mastectomy  specim  en  r  esults  at  higher  dose  rev  ealed  that  for  both  MCI  and  tomosynthesis, 
highest  dose  setting  and  m  aximum  angular  span  with  an  angular  separation  of  2.75  0  m  ay  be 

optimum,  indicating  a  threshold  in  the  num  ber  of  projections  per  angular  span  for  optim  um 
performance.  Most  importantly,  the  study  demonstrated  that  the  performance  of  optimized  breast 
MCI  may  exceed  that  of  mammography  and  tomosynthesis  by  18%  and  8%,  respectively. 
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OVERALL  CONCLUSIONS 

In  this  study,  an  optimization  firam  ework  was  developed  to  m  aximize  the  diagnostic 
performance  of  MCI.  It  was  revealed  that  the  peak  performance  for  MCI  at  the  clin  ically  relevant 
dose  levels  of  one-  and  two-view  mammography  may  be  achieved  at  15  -  17  projections  spanning 
an  angular  arc  of  ~45  °,  the  widest  angle  tested  in  this  st  udy.  Most  importantly,  we  dem  onstrated 
feasibility  ofm  ulti-plane  cor  relation  im  aging  ( MCI)  as  a  technique  f  or  im  proved  breast  can  cer 
detection.  Compared  to  mammography  and  tomosynthesis,  MCI  was  shown  to  be  potentially  m  ore 
accurate,  and  cost-  and  dose-  effective.  T  his  study  thus  paves  way  for  future  clinical 
implementation  of  MCI. 
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